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Abstract. Particle acceleration is expected to take place at shocks that form during 
the process of large scale structure formation. Electrons accelerated at such shocks 
can upscatter a small fraction of the photons in the cosmic microwave background 
up to the gamma ray band. Here we make predictions about the detectability of 
the 7-ray emission from forming clusters of galaxies with future GeV and TeV 
gamma ray telescopes. We also estimate the contribution of these sources to the 
extragalactic diffuse gamma ray background. 



1. Introduction 

Non-thermal emission from ric h clus- 
ters is observed in the radio llFerettiJ 
Ithese proceedings and references therein*) 
and hard X-ray ^usco Femiano ct al. 
Il999j) bands while only upper limits for 
their gamma ray fluxes have been obtained 
from an accurate an alysis of EGRET data 
ijReimer et a1. 11200^. 

Despite the fact that to date none of 
the nearby rich clusters has been detected, 
there are theoretical arguments suggest- 
ing that clusters could be sources of high 
ener gy radiation. iBerezinskv et aL I l|l997j) 
and IVoelk et al. I ^^199(1^ first understood 
that cosmic ray protons accelerated in the 
intracluster medium (ICM) would be con- 
fined there for cosmological time scales, en- 
hancing the possibility of inelastic proton- 
proton collisions and consequent gamma 
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ray production through neutral pion de- 
cay. Moreover, shocks associated with large 
scale structure formation might accelerate 
electrons up to TeV energies. These elec- 
trons can upscatter a small fraction of the 
photons in the cosmic microwave back- 
ground radiation up to gamma ray energies 
through inv erse Compton scattering (ICS) 
ijLoeb fc Waxman 200Q) . For these reasons 
the detection of gamma ray emission from 
the ICM represents one of the major goals 
of future gamma ray telescopes. 

Here we addressed the issue of the de- 
tectability of clusters in gamma rays and of 
the contribution of large scale structure for- 
mation to the extragalactic diffuse gamma 
ray background (EDGRB). We restrict our 
attention only to ICS of rclativistic elec- 
trons accelerated at shocks generated dur- 
ing structure formation. The contribution 
from protons to the gamma ray flux can 
be of great importance, but it is also quite 
uncertain because, due to confinement, the 
amount of protons diffusively stored in the 
ICM depends on the whole history of the 
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cosmic ray sources in clusters. Any con- 
tribution to the gamma ray emission from 
protons can only increase the fluxes derived 
below, that can be considered as conserva- 
tive estimates. 

All the presented results are obtained 
by means of semi-analytical calculations, 
while an approach based on numerical sim- 
ulations is described by Miniati (these pro- 
ceedings). 

2. Shock waves during structure 
formation: mergers and accretion 

Clusters of galaxies form hierarchically 
through the merger of smaller dark mat- 
ter halos. During mergers shock waves nat- 
urally form and heat the ICM up to the 
observed temperatures. 

A useful analytical description 
of this process was proposed by 
iPress fc Schccht cr (1974) and later devel- 
oped bv iLacev fc Cole I |i99^ in the form 
of the so called extended Press-Schechter 
formalism. This approach provides us with 
the comoving number density of clusters 
with mass M at cosmic time t, n{M, t) 
and with the rate at which clusters of mass 
M merge at a given time to form clusters 
with final mass M', n{M,M',t). It is 
easy to use this formalism to construct 
simulated merger trees for clusters with a 
given present mass. 

Tn ICabici Riasi I l|2nn8all we devel- 
oped an analytical recipe that allows us to 
estimate the Mach number of merger re- 
lated shocks, once the masses of the two 
merging clusters and the redshift at which 
the merger event occurs are given. The rel- 
ative velocity Vr can be evaluated from en- 
ergy conservation: 

_ GM,M2 ^ l^^y ^ GA'hM2 
Ri + R2 2 2Rta 

where Mi and Ri are the masses and the 
virial radii of the two merging clusters, 
Mr is the reduced mass and Rta is the 
turnaround radius, where the two clusters 
are supposed to be at rest. The sound speed 
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Fig. 1. Cumulative distribution of Mach 
numbers in mergers o f clusters of galaxies 
llCabici & Blasi ll2004D . 

in the z-th cluster, needed in order to eval- 
uate the shock Mach number, follows from 
the virial theorem: 

cl, ^ GMJ{2R,) . (2) 

The procedure illustrated above can be 
applied to a generic merger event in the 
simulated history of a cluster. The results 
obtained for 500 different realizations of the 
merger tree for a cluster with present mass 
IO-'^^Mq are plotted in fig. ^ where the cu- 
mulative probability distribution of shock 
Mach numbers is shown. 

The striking feature of this plot is the 
fact that most merger related shocks are 
mildly supersonic, having Mach numbers ~ 
1.5, while only 6% of the s hocks have Mach 
numbe r greater than 3 l|Gabici &: Blasi I 
120081. 

Parallel to the merger processes, a sec- 
ondary infall of gas occurs at all times 
onto the potential well which is being 
formed. The information about the virial- 
ization of the inner region of this accre- 
tion flow is propagated outwards through 
a so called accretion shock. While merger 
related shocks form and propagate in the 
hot virialized ICM, accretion shocks propa- 
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gate in the cold (W^-IO^K), non virialized 
intergalactic medium, an d therefore they 
are exp ected to be strong l|Gabici fc Blasi I 

Mm. 

There is still some debate on the typi- 
cal Mach number of the shocks developed in 
the ICM during cluster mergers. In fact, re- 
sults obtained by means of numerical simu- 
lations suggest the presence of a large frac- 
tion of hi gh Mach number sho cks in cluster 
mergers llMiniati et al. I| 2000l). More recent 
simulations (iT^v u et a,l. 1 12003'^ seem to find 
more weak shocks than in Miniati et al. I 
l)2000j) . However the comparison with the 
analytical calculations presented here ap- 
pears difficult because of a different classifi- 
cation of the shocks in the two approaches. 
On the other hand, the weakness of the 
merger-related shocks seems to be also sug- 
gested by the few observations in which the 
Mach number of the shock can be measured 
fe.g..[MaTkevitch et aP |2003)). 

The strength of the shocks is of cru- 
cial importance for the acceleration of 
suprathermal particles and related emis- 
sion, as discussed in the next section. 



3. Shock acceleration during 
structure formation 

Diffusive acceleration can take place at 
merger and accretion shocks and a small 
fraction of the particles in the accreting gas 
is estracted from the thermal distribution 
and energized up to relativistic energies. 

The maximum energy achievable for a 
particle can be estimated equating the ac- 
celeration time: 



AD{E) 



(3) 



where D{E) and u are the energy 
dependent diffusion coefficient and the 
flow spee d, b oth e valuated upstream 
l)Gabici & Blasi 1l2004l) . with the minimum 
between the energy loss time, the escape 
time or the life time of the accelerator. 
Since we are considering electrons, the rel- 
evant time scale is the energy loss time due 



to ICS off the microwave background radi- 
ation photons: 



E 



Tics 



^arcucmbl^ 



(4) 



where Ucmb is the energy density in the 
CMB radiation and the other symbols have 
their usual meanings. The highest value of 
the maximum energy is obtained adopting 
the Bohm diffusion coefficient: 



Er, 
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lO^cm/s 



TeV .(5) 



At these energies the ICS off the CMB pho- 
tons occurs in the Thomson regime, there- 
fore the maximum energy of the radiated 
photons is given by: 



Ery 



(6) 
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which falls in the very high energy gamma 
ray band. It is important to stress, how- 
ever, that the Bohm diffusion coefficient 
is the smallest possible one, and that a 
different choice of the diffusion coefficient 
could result in a too low maximum en- 
ergy of the electrons, so that the electron 
population would be unable to upscatter 
the CMB photo ns up to gamma ray en- 
ergi es.Following iLoeb. &:_Waxman (200Qi) 
an d iGabici fc Bl asi ( 20033), we will use for 
the maximum energy the value obtained in 
eq. El 

Particle spectra can be easily calculated 
in the test-particle limit, when the back re- 
action of the accelerated particles on the 
shock structure can be neglected. In this 
case the particle distribution function at 
suprathermal energies is a power law in mo- 
mentum N(j>) oc with slope connected 
to the shock Mach number A4 through the 
well known relation: 



a 



(7) 
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It follows that, in order to accelerate 
flat particle spectra and to have apprecia- 
ble gamma ray emission strong shocks are 
needed. As seen in the previous section, ac- 
cretion shocks are strong, while the major- 
ity of merger related shocks are weak and 
accelerate very steep particle spectra, irrel- 
evant from the point of view of high energy 
phenomena related to electrons. 

4. Clusters as 7-ray sources and 
their contribution to the EDGRB 

The number of accreting clusters detectable 
by a gamma ray telescope with sensitivity 

^lirn is. 

NacciFum) - (8) 

r°° dv r°° 

/ dz— / dMn{M,z) 

Jo dz JM(Fii^,z) 

where V is the comoving volume between 
redshift z and z-\-dz, n{M, z) is the comov- 
ing density of objects as a function of clus- 
ter mass and redshift. The mass of a cluster 
accreting at redshift z whose flux is Fum, 
M{Fiim, z), is evaluated assuming that the 
accretion shock is located at the cluster 
virial radius and that 5% of the kinetic en- 
ergy crossing the shock sur face is converted 
into relativistic electrons (jGabici fc Blasi I 
|2003b, 2004). 

For merging clusters a similar ex- 
pression, involving the merger rate 
Tl{M,M',z), holds. The only difference is 
that for accreting clusters the spectra of 
the accelerated particles are power laws 
with fixed slope a = 2 while for merging 
clusters the spectral indexes are evaluated 
self-consistently as described in ^ 

The number of clusters with gamma ray 
emission at energies in excess of lOOAfeV^ 
larger that some fixed flux F is plotted in 
fig- HI The solid and dashed lines refer to 
accreting and merging clusters respectively. 
The sensitivity of the EGRET, AGILE 
and GLAST instruments are also shown. 
According to our calculations GLAST and 
AGILE should detect a few tens and ^ 
10 clusters respectively. Our results are 
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Fig. 2. Number of accreting (solid line) 
and merging (dashed line) clusters with 7- 
ray flux greater than F as estimated by 
iGab ici & Blasi (20 (W). Open circles refe r 
to the simulation bv lKeshet et al. I l)2003j) . 



in agreement with the lack of associa- 
tion between cl usters and EGRET sources 
l|R,eimer et al. l l2003^. 

The superposition of the emission from 
single clu sters results in a diffuse back- 
ground. In lGabici fc Blasi ( 2003b) we esti- 
mated the contribution to this background 
from forming structure and we summarize 
our results in figEJ The shaded yellow area 
represent s the EDGRB as me asured by 
EGRET l|Sreekumar et al. Ill998t) while the 
green and red lines represent the contribu- 
tion to the background from merging and 
accreting clusters. A more realistic estimate 
for the contribution of merging clusters can 
be obtained assuming a minimum mass for 
the merging halos of IO^^Mq (blue line). 

Since our results greatly differ 
from some previous results, some 
comments and compariso ns ar e re- 
quired. IWaxm an & Locb l|2000D and 
iTotani fc Kitavama I (\2Q()(h developed 
a Press-Schechter based method similar 
to the one presented here. The authors 
claimed that at least a few tens of clusters 
should have been visible to EGRET. 
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Fig. 3. Diffuse gamma ra y emission from 
structu re formation l|Gabici fc Blasi I 
l)2003bD . see text for details). 



GLAST, on the other hand, was predicted 
to be able to detect more than a few 
thousands of such objects. This predic- 
tion seemed to be supported by some 
preliminary observational evidences for an 
association between unidentified high lati- 
tude EGRET sources and A bell clusters (or 
cluster pairs ) presented bv iGolafrancesco I 
l|2002|) and iKawasakv fc Totani I l)2000(l . 
The statistical significance and physical 
plausibility of such an as sociation was 
strong ly questioned by iReimer et al. I 
l)2003() . in which the important conclusion 
was reached that we still have to await 
the first observational evidence for the 
high-energy gamma-ray emission of galaxy 
clusters. These findings seem to be per- 
fectly in line with the predictions described 
in the present paper. 

All th e shocks were assumed t o be 
strong bv iTotani fc Kitayama I l|2000D and 
IWaxman fc Loebl l)20fl0t ). As a conse- 
quence, the spectrum of the accelerated 
particles was always taken to be oc p^^. 
As discussed ab ove, and as d escribed in 
great er details bv lGabici fc Bla si ( 2003^3, 
l2nnl . for merger related shocks this as- 
sumption leads to incorrect results. The 



gamma ray emission is overestimated by or- 
ders of magnitude and its spectrum does 
not reflect the real strength of the shocks 
developed during mergers of clusters of 
galaxies. Similar arguments can be used 
in order to explain the difference between 
our estimate of the dif fuse emis sion and 
the results presented bvlLoeb fc, W axmani 
^ l|2000(l and lTotani fc Kitavama n 2000). in 
which the whole EDGRB was saturated by 
the contribution from structure formation. 

Our findings are in good agreement 
with the results from the num erical simula- 
tions by Keshe t et al. 1 1|2003() (open circles 
in fig. 12 and in which only gamma rays 
from ICS were considered and the same 
assumption we made on acceleration effi- 
ciency was used. 



5. TeV emission from clusters 

The maximum energy of the upscattered 
CMB photons falls in the TeV range (see 
eq|n|) so that we can wonder if clusters could 
be detected by ground based Cherenkov 
telescopes. In fig. ^ we plot the sensitiv- 
ities for a generic telescope as calculated 
in Aharonian ct al. (1993) for 100 hours 
of observation with an array consisting of 
10 cells (pointlike and 1° sources) together 
with the expected spectra for three differ- 
ent cases (top to bottom): 1) a merger be- 
tween two clusters with masses 10^^ and 
IO^^Mq, 2) a IQI^Mq accreting cluster 
with a magnetic field in the upstream re- 
gion equal to 0.1/xG, 3) the same as 2) 
but with a ten times lower magnetic field. 
All the considered clusters are at a dis- 
tance of lOOMpc. The thin continuous and 
dotted lines refer to the spectra calcu- 
lated with and without taking into account 
the absorption on the infra red background 
(|De Jager fc Stecker ll200^ . 

The apparent size of a IO^^Mq accret- 
ing cluster at a distance of lOOMpc is com- 
parable or even greater than the field of 
view of Cherenkov telescopes a few de- 
grees) , while the size of the emitting region 
for the merger considered in fig. 01 should 
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Fig. 4. Gamma ray emission in the TeV 
region. The thick sohd hnes represent the 
sensitivities of a cherenkov array for point 
sources (lower curve) and extended sources 
(upper curve). The predicted gamma ray 
fluxes from a Coma-like cluster with and 
without absorption of the IR-background 
are plotted as dashed and solid lines. 

be roughly 1°, making it marginally de- 
tectable. 

6. Conclusions 

We developed a recipe to evaluate the Mach 
numbers of th e shocks that form dur ing a 
merger event ijGabici fc Blasi ll2003a^ . We 
found that these shocks are often weak 
and the corresponding spectra of shock 
accelerated particles ar e ste ep (see also 
iBerrington fc Dermer I ll2003tl). We gen- 
eralized our approach in iGabici fc Blasil 
(|20D3b','200?) where we considered also the 
strong accretion shocks and we made pre- 
dictions about the number of clusters de- 
tectable in gamma rays by future gamma 
ray telescopes and about their contribution 
to the EDGRB. 

According to our prediction, AGILE 
and GLAST should be able to detect ~ 10 
and a few tens of objects while the dif- 
fuse background produced by clusters is ex- 



pected to be '--^ 10 times smaller than that 
observed by EGRET. 

The detection of clusters in the TeV 
energy range could be problematic due to 
their extended apparent size. 
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